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Abstract 

System analysis of oxygen production by hydrogen reduction of lunar regolith has shown the 
importance of the relative time scales for regolith heating and chemical reaction to overall performance. 
These values determine the sizing and power requirements of the system and also impact the number and 
operational phasing of reaction chambers. In this paper, a Nusselt number correlation analysis is 
performed to determine the heat transfer rates and regolith heat up times in a fluidized bed reactor heated 
by a central heating element (e.g., a resistively heated rod, or a solar concentrator heat pipe). A coupled 
chemical and transport model has also been developed for the chemical reduction of regolith by a 
continuous flow of hydrogen. The regolith conversion occurs on the surfaces of and within the regolith 
particles. Several important quantities are identified as a result of the above analyses. Reactor scale 
parameters include the void fraction (i.e., the fraction of the reactor volume not occupied by the regolith 
particles) and the residence time of hydrogen in the reactor. Particle scale quantities include the particle 
Reynolds number, the Archimedes number, and the time needed for hydrogen to diffuse into the pores of 
the regolith particles. The analysis is used to determine the heat up and reaction times and its application 
to NASA’s oxygen production system modeling tool is noted. 


Nomenclature 

A reactor cross section area (m 2 ) 

A e f f effective surface area of heater (m 2 ) 

Ar Archimedes number 

c molar concentration of hydrogen in particle (moles/m 3 ) 

c s molar concentration of hydrogen at reaction surface (moles/m 3 ) 

c x molar concentration of hydrogen at location x in reactor (moles/m 3 ) 

c 0 molar concentration of hydrogen at upstream end of regolith bed (moles/m 3 ) 

C e ff effective specific heat (J/kg/K) 

C pg specific heat of gas (J/kg/K) 

C ps specific heat of regolith (J/kg/K) 

D effective gas diffusion coefficient (m 2 /s) 

E energy required to heat bed 

h heat transfer coefficient 

g gravitational acceleration (m/s 2 ) 

k equilibrium constant 

L length of regolith bed (m) 

M c it effective mass (kg) 
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A/t 0 tai total molar water yield (mol) 

N number of regolith particles in the bed 

Nu Nusselt number 

Pr Prandtl number 

r p regolith particle radius (m) 

Re Reynolds number 

t time (s) 

t c chemical time (s) 

t D intra-particle diffusion time (s) 

tn l imi t chemical time in diffusion limit (s) 

t R reactor residence time (s) 

tftasLimit chemical time in equilibrium limit (s) 

T temperature (K) 

T i, heater temperature (K) 

T u T 2 initial and final temperature of bed (K) 
u velocity through reactor (m/s) 

Rbed volume of regolith bed ( m 3 ) 
x axial coordinate (m) 

a particle conversion 

ao particle conversion at x = 0 

P regolith thermal diffusivity (m 2 /s) 

8 thermal conduction length scale (m) 

s void fraction in regolith bed 

p a molar concentration of ilmenite in regolith (mol/m 3 ) 

p„ gas density (kg/m 3 ) 

p. s regolith density (kg/m 3 ) 

p gas viscosity (Ns/m 2 ) 

X gas thermal conductivity 

¥ molar outflow of water (mol/s) 


Introduction 

Hydrogen reduction of lunar regolith has been proposed as a viable technology for oxygen production 
on the Moon. (Taylor and Carrier 1992, McKay and Allen 1996, Rice, Hermes and Musbah 1997). The 
reduction processes involve solid-gas reactions. Such systems are well developed for earth applications 
but mainly for the purpose of metal extraction as opposed to water vapor and oxygen production for lunar 
applications. In addition, in the lunar context, the reactor must interface with the other sub-system 
processes such as upstream regolith extraction and beneficiation, and downstream electrolysis and phase 
separation in a way that establishes the most favorable balance between efficiency, robustness, and 
equivalent system mass. For example, the reduction of ilmenite may be represented as: 

FeO.Ti0 2 + H 2 <=> Fe + Ti0 2 + H 2 0 (1) 

It is carried out in the reactor sub-system. The formed water may be electrolyzed in the electrolyzer 
subsystem to produce oxygen and hydrogen. 

From upmass and cost considerations, it is essential to operate these systems on the lunar surface in as 
nearly a closed-loop manner as possible for the hydrogen. It is in this context that knowledge of 
processing times in the various sub-systems becomes important. For example, a mismatch between the 
flux of water exiting the reactor and the amount of water that may be processed by the electrolyzer must 
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be mitigated by design such as sizing the appropriate buffer tank. This paper investigates the factors that 
influence the overall processing time in the reactor sub-system. 

The temperature of the reaction is one of the important parameters of the process, llmenite is just one 
constituent of mare basalt and it cannot be assumed that regolith beneficiation (e.g., in this case, removal 
of the nonilmenite components) will be carried out prior to the reactor processes. The sintering and 
melting temperature of the unbeneficiated regolith is noted to be in the range of 1370 to 1670 K 
(Chambers 1995, McKay and Allen 1996). However, hydrogen reduction has been found to be effective 
also at lower temperatures, e.g., in the range of 970 to 1270 K (Yoshida et al., 2000). Therefore, while 
three-phase (gas-liquid-solid) reactors for hydrogen reduction may be applicable at post-melt 
temperatures, these are not considered here and attention is restricted to gas-solid interactions. In order to 
provide a basis for discussion, it is assumed that the chemical reduction of the regolith occurs in a gas- 
solid fluidized bed reactor. 

The conversion process in the reactor consists of the following overall steps (1) transfer of regolith 
into the reactor, (2) fluidization and heating of the regolith, (3) reaction of hydrogen with the ilmenite 
portion of the regolith, (4) outflow of left-over hydrogen and produced water vapor, and (5) removal of 
the spent regolith from the reactor. In this paper, attention is focused on steps (2) to (4). 

The reported values for the thermal diffusivity, |3, of lunar regolith are very low (Heiken, Vaniman 
and French 1991). In the absence of a gas flow, the time to heat the regolith bed over a distance 8 is 
proportional the ratio 8 2 /p, which results in unacceptably large values. On the other hand, hydrogen has a 
relatively high thermal conductivity and helps to stir the bed so that having hydrogen flow through the 
bed decreases the heating time appreciably. Thus, it is advantageous to fluidize the bed during the regolith 
heating step. There are also advantages related to reducing regolith particle agglomeration and sintering 
which can adversely impact reaction times and water production rates. 

The objectives of this paper are (1) to develop analytical models for regolith heating and for chemical 
reduction by hydrogen in order to identify parameters influencing the heating and reaction times, (2) to 
identify the parameters that govern these times, and (3) to obtain relationships for the reaction time in the 
two limits of diffusion (i.e., mass transfer) control and chemical equilibrium control. As noted earlier, 
these times are important from the overall oxygen production system perspective. The thermal analysis is 
based upon a Nusselt number correlation. The reaction model considered is based upon a “shrinking core” 
diffusion-reaction model for the gas-solid reduction process (Levenspiel 1998) which is coupled with the 
inflow/outflow of gases. The methodology has been discussed in our earlier work (Hegde, 
Balasubramaniam and Gokoglu 2007). It appears consistent with results of regolith particulate analysis 
conducted on Apollo-returned samples (McKay and Allen 1996) which suggest that the ilmenite particles 
are not completely solid but have some porosity that would enable gases to diffuse into the particles, i.e., 
the FeO portion can be essentially completely converted leaving behind the reduced iron. 

In the next section, the regolith heating is considered first and a model for the heating time and 
required power is developed. Pertinent results from the shrinking core analysis are also presented. 
Subsequently, in the Results section, factors impacting these processes are evaluated and discussed and 
applications of the model to NASA’s oxygen production system modeling tool are noted. 

Model Development 

In the first sub-section below, the time taken to heat the regolith is calculated using a Nusselt number 
correlation for fluidized bed heat transfer. This is followed by a summary of the particle conversion 
analysis and the coupling between the particle conversion and the bulk reactor flow. 

Heat Transfer Analysis 

It is assumed that the regolith is heated by a cylindrical rod located centrally in the cylindrical reactor. 
The surface temperature of the heater rod can be controlled e.g., by adjusting the current through it as in 
an electrical heater or by controlling the energy flux as in a solar concentrator heat pipe. The different 
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pathways for heat transfer in the bed are gas convection and conduction, particle motions, and possibly, 
radiative transfer. 

Fluidized bed heat transfer correlations developed in the literature under convective-conductive 
conditions suggest that the controlling step is the heat transfer from the heater surface to the gas film next 
to it. Direct transfer of heat to the bed particles adjoining the heater is not efficient because the contact 
area between particles and the heater surface is small. If the bed is well stirred by fluidization, its 
temperature is nearly uniform. In such a case, a Nusselt number may be defined which relates the heat 
transfer to the bed through the gas film in contact with the heater, i.e., the regolith bed is taken to be a 
uniform entity. The reactor wall is not usually insulated from the bed, so in general, there will also be heat 
transfer from the bed to the reactor wall. If the reactor wall is not too thick and its thermal conductivity is 
high, the wall may be assumed to be at a uniform temperature. Furthermore, if the reactor’s thermal mass 
(i.e., the product of its mass and specific heat) is small compared to that of the bed, it is reasonable to 
assume that the temperature of the reactor is the same as that of the bed. 

The heating rate with the above scenario is given by: 

M etf C eff ~7T = { T h ~ T ) ’ ( 2 ) 

at 

where the thin film heat transfer coefficient, h, is related to the Nusselt number, Nu, by: 

Nu = 2hr p / X . (3) 


In Equation (2), M e{f C t ff is the effective value of the product of the mass and specific heat and is a 
weighted average of the values for the gas, regolith, and reactor. A eff is the effective or active surface area 
of the heater rod in contact with the bed. 

The solution of Equation (2) for the instantaneous bed temperature T is given by: 

T = T h - (T h - 7j ) exp (- hA eS t /(m eff C eff )) . (4) 

Conversely, the time to raise the bed temperature from 7), to a temperature T 2 is given by: 


{ h ~ 


^effQsff | n T h ~ T 2 
hd c ff ^ T h - 7j y 


(5) 


Various fluidized bed heat transfer Nusselt number correlations are available in the literature. For 
example, the following one (Borodulya et al. 1991) appears to describe the situation adequately: 


Nu = 0.74Ar 


o.i 


f A 

P, 

0.14 

f c ) 

^ ps 



r 

Vps) 


0.24 


(l - s) 2 ' 3 + 0.46 Re Pr 


Mi! 


( 6 ) 


The Archimedes number, Ar, the Reynolds number, Re, and the Prandtl number, Pr, are defined as 
follows: 


Ar = 8g7 ,3 pg(p ? - Pg)/p 2 , 


(7) 


Re = 2ur p p„ / p 


( 8 ) 
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¥x = C pg \i!'k g . 


( 9 ) 


The energy, E, required to heat the bed and reactor from T\ to T 2 , in the absence of other losses (e.g., 
radiative), is given by: 


E = M ef{ C ef{ (T 2 -T l ). (10) 

An average power required is obtained by dividing E by the time, T h . Note that the instantaneous power, 
M e ffC e ffDT/dt, can be obtained by substituting the solution for T (Eq. (4)) into Equation (2). 

Since radiative losses vary as the fourth power of temperature, they can become important as the 
reactor is heated to higher temperatures. It is expected that there would be insulation around the reactor 
that would reduce this loss. This loss can be accounted for by including it in the energy required. 

Chemical Conversion Analysis 

The chemical conversion analysis couples the “shrinking core” model with the inflow and outflow of 
the gases and has been discussed in earlier work (Elegde, Balasubramaniam and Gokoglu, 2007). Elere, 
some results are presented to highlight the findings from the model and to show how the characteristic 
times of the process are obtained. 

Defining the conversion a by: 


a = 1 



( 11 ) 


the shrinking core analysis of a single particle yields a relationship between a and the associated 
conversion time, t\ 


3[l - a] 2 ' 3 - 2[l - a]- 1 = _ . 

P a r p 


( 12 ) 


Further, assuming equilibrium conditions at the reaction surface, the above equation reduces to: 

3[i-„p-2[i-q]-l=- 60fa A. (13) 

d + *)p „l> 2 

Note that an expression for the equilibrium constant, k, as a function of only the reaction temperature has 
been provided (William and Erstfeld, 1979): 

k _ io(- 2126 - 1/r+0 - 6439 ) _ (14 ) 

The time for complete conversion is given by setting a to unity. 
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Coupling Particle Analysis With the Reactor Bulk Flow 

Assuming that the gas velocity, concentration of hydrogen and the associated conversion a are 
functions of axial distance along the reactor, and the inlet hydrogen concentration is constant in time, the 
following equation has been obtained: 


-4 In 


1 — (l — a) 1/3 
l-(l-a 0 ) 1/3 


+ In 


2(l-a) 1/3 +l 

2(l-a 0 ) 1/3 +l 



(15) 


where ao is the value of a at x = 0 and is given by Equation (13). 

In certain cases, closed form solutions may be obtained for the concentration of hydrogen and extent 
of conversion in the reactor. An example is when the regolith and gases are sufficiently stirred or well- 
mixed so that gradients in conversion and gas concentration are smeared out. In this well-stirred limit, the 
concentration of hydrogen exiting the reactor is also the average concentration in the reactor and is given 

by: 


'M= 


c 0 

1 + k 


1+1 


1 + 3(1 - s) 


f \ 

D 

ji 


1-a l/ ' 
1 - (l - a) 1 


(16) 


Results 

Regolith batch sizing will depend upon the oxygen production required and the processing time per 
batch. In general, the processing time per batch is the sum of the heating and chemical conversion times. 
This section considers factors influencing these times based upon the analysis presented earlier. 

Heating Time 

Equation (5) indicates that the heating time is related directly to the Nusselt number and batch size 
(i.e., regolith mass) and inversely related to the active surface area of the heater. The Nusselt number 
depends upon the Reynolds number of the gas flow relative to the regolith particles, the Archimedes 
number, and the Prandtl number of the fluidizing gas. The Prandtl number is practically constant with 
temperature. The Reynolds number is related to the fluidizing velocity and particle size. The particle size 
is fixed by the regolith characteristics and there may be some flexibility in this regard depending upon the 
preprocessing facilities available. The fluidizing velocity is also related to the particle size. The 
Archimedes number is mainly dependent upon the regolith density (which is much greater than the gas 
density). The most flexible way of changing the heating time seems to be by changing the regolith batch 
size, active heater area, and/or power input to the heater. 

The active heater surface area may be increased by increasing the heater diameter. For a fixed reactor 
diameter, this also increases the height, L, of the regolith bed which is also conducive to heat transfer. 
However, this may have implications for the reaction time as will be discussed later. The regolith batch 
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size may be reduced to decrease the heating time but this may be restricted from overall oxygen 
production considerations. 

The power input to the heater is related to the heater temperature T h . From Equation (5), it is clear that 
increasing the heater temperature will decrease the time to heat the regolith to a given temperature T 2 . 
Flowever, there are limitations to raising 7), substantially over T 2 e.g., the presence of hot spots or regions 
in the regolith bed which are at higher than desired temperature. 


The shrinking core analysis provides an expression for the conversion time of a single particle. In 
many instances it is assumed that the conversion time of the regolith bed is the same as the single particle 
conversion time. Equation (15) provides a means for determining the conditions under which this holds. 
Choosing x = L (i.e., the height of the regolith bed), it is clear that the right hand side of Equation (15) is 
related to the ratio of the gas residence time in the bed, Liu, to its diffusion time in a particle, r p 2 /D. This 
time factor ratio plays an important role in determining the extent of conversion in the bed. For example, 
the upstream end (entrance region) of the bed is always exposed to fresh hydrogen and would be the first 
to get converted. The time ratio factor gives a measure of the conversion in the remaining portions of the 
bed. Figure 1 plots the conversion, a (L), at the downstream end (exit region) as a function of the time 
ratio factor for the case where ao is equal to unity, i.e., when the particles at the entrance are already 
completely converted. As the time ratio factor becomes small (i.e., the intra-particle diffusion time is 
much greater than the gas residence time in the regolith bed) it can be seen from Figure 1 that the 
conversion, a (L), also approaches unity, i.e., the particles at the downstream end of the bed (and more 
generally, at any location in the bed) are also consumed fully. This indicates that the time rate of chemical 
conversion for all particles in the bed is given by Equation (13). The corresponding time for complete 
conversion may be termed the chemical time in the diffusion controlled limit. Setting a to one in 
Equation (13) yields: 


Chemical Time 



6 DkcQ 


(17) 


25 


20 - 


\ 


15 - 


Time Ratio Factor, 10 


3(1 -£)DL/(r p 2 u) 


5 - 


\ 


0 
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0.1 


0.2 0.3 0.4 0.5 0.6 0.7 


0.8 


0.9 


a (L) 


Figure 1. — Dependence of the Regolith Conversion at the Downstream End of the Reactor 
Bed Upon the Time Ratio Factor 3(1 -s) DL/r p 2 u. 
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On the other hand, as the residence time becomes much larger compared to the diffusion time, i.e., the 
time ratio factor becomes large, it can be seen from Figure 1 that the conversion occurs first primarily 
near the upstream end of the bed and not much conversion occurs downstream since a (L) is small for this 
condition. This is referred to as the chemical equilibrium limit, i.e., the residence time in the bed is long 
enough so that the water concentration reaches its chemical equilibrium value in the gas phase. At this 
point no further conversion takes place. In this limiting case, the chemical time may be estimated by 
assuming equilibrium concentration of water in the outflow. Then, the molar outflow of water is given by: 

£ 

^ = Y+k) UACo '‘ (18) 


and the total molar water yield is: 


^ total 'I^ResLimit • (19) 

Equating this to the available molar water yield in the regolith provides the chemical time in the residence 
time limit or equilibrium limit: 


^ResLimit Pa(l e )^bed/V- (20) 

Since the conversion time in the diffusion limit depends inherently upon the particle characteristics such 
as its size and the intra-particle gas diffusivity, it is clear that this is a lower limit to the conversion time. 
The equilibrium limiting conversion time is variable however, depending upon the flow velocity and bed 
height. Varying these quantities is a way to slow down the conversion rate. This may be useful from an 
overall system design perspective. 


Application to System Design 

The analysis of heating and chemical times presented here may be applied to sizing and design of the 
reactor assembly. It has been noted in earlier studies that the regolith reduction rate increases directly with 
the temperature, (McKay and Allen 1996) presumably due to increase in the effective diffusion 
coefficient (Hegde, Balasubramaniam and Gokoglu 2007). This implies that most of the conversion 
occurs at the operating temperature. In other words, the heating and chemical reductions are essentially 
serial operations for a given batch of regolith. Thus, a system comprising two reactors may be used for 
increased throughput. Ignoring for current purposes, the regolith loading and unloading times into the 
reactor chamber, it is clear that designing for equal heating and reaction times will ensure continuous 
production of oxygen, i.e., as one batch of regolith is being heated the other batch is being chemically 
reduced. 

The equations provided here provide a framework for achieving the proper design parameters. The 
diffusion limited chemical time (Eq. ( 1 7)) provides the minimum time for conversion. Equations (5) and 
(10) in conjunction provide the power and heater rod characteristics (e.g., effective area) for a given 
heating time, e.g., the minimum reaction time. If the power available is less than the heating power 
required, the heating time has to be increased. In such a case, the reaction time needs to be modified e.g., 
by driving the reactor dimensions toward the residence limited conversion time (Eq. (20)). Clearly, this is 
an iterative process. Currently, NASA is developing an overall system model incorporating these ideas 
(Santiago-Maldonado 2007). 
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Conclusions 


Analytical models for regolith heating and for chemical reduction by hydrogen have been developed 
to identify parameters influencing the heating and reaction times. These times are important from the 
overall oxygen production system perspective. The regolith heating model is based upon a Nusselt 
number correlation for fluidized bed heat transfer. While particle size and parameters such as the particle 
Reynolds and Archimedes numbers play a role in the heating time determination, it appears that the most 
feasible way of controlling this time is by adjusting the active heater surface area and mass of the regolith 
per batch. The chemical reduction model results from coupling a shrinking core type of analysis with the 
overall regolith bed flow characteristics. It is shown that the chemical time in the diffusion limit is also 
the lower limit for the conversion time. However, by adjusting the residence time of the gases in the 
regolith bed, the conversion time can be increased. 
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